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Abstract 



Gamma correction can be expressed by the equation 



This paper presents invariants under gamma correction and 
similarity transformations. The invariants are local features 
based on differentials which are implemented using deriva- 
tives of the Gaussian. The use of the proposed invariant rep- 
resentation is shown to yield improved correlation results in 
a template matching scenario. 



1 Introduction 

Invariants are a popular concept in object recognition and 
image retrieval [jjj, ||, [7|, [14], |T5j ], They aim to provide 
descriptions that remain constant under certain geometric or 
radiometric transformations of the scene, thereby reducing 
the search space. They can be classified into global invari- 
ants, typically based either on a set of key points or on mo- 
ments, and local invariants, typically based on derivatives of 
the image function which is assumed to be continuous and 
differentiable. 

The geometric transformations of interest often include 
translation, rotation, and scaling, summarily referred to as 
similarity transformations. In a previous paper Jl2|], build- 
ing on work done by Schmid and Mohr [ 1 1 1, we have pro- 
posed differential invariants for those similarity transforma- 
tions, plus linear brightness change. Here, we are looking at 
a non-linear brightness change known as gamma correction. 

Gamma correction is a non-linear quantization of the 
brightness measurements performed by many cameras dur- 
ing the image formation process.[] The idea is to achieve bet- 
ter perceptual results by maintaining an approximately con- 
stant ratio between adjacent brightness levels, placing the 
quantization levels apart by the just noticeable difference. 
Incidentally, this non-linear quantization also precompen- 
sates for the non-linear mapping from voltage to brightness 
in electronic display devices . 



(1) 



Historically, the parameter gamma was introduced to describe the non- 
linearity of photographic film. Today, its main use is to improve the output 
of cathode ray tube based monitors, but the gamma correction in display 
devices is of no concern to us here. 



where / is the input intensity, I 1 is the output intensity, and 
p is a normalization factor which is determined by the value 
of 7. For output devices, the NTSC standard specifies 7 = 
2.22. For input devices like cameras, the parameter value is 
just inversed, resulting in a typical value of 7 = 1/2.22 = 
0.45. The camera we used, the Sony 3 CCD color camera 
DXC 950, exhibited 7 w 0.6]] Fig. [j] shows the intensity 
mapping of 8-bit data for different values of 7. 




Figure 1: Gamma correction as a function of intensity, 
(solid) 7=1; (dashed) 7 = 0.45; (dotted) 7 = 2.22. 
Note how, for 7 < 1, the lower intensities are mapped onto 
a larger range. 



It turns out that an invariant under gamma correction can 
be designed from first and second order derivatives. Ad- 
ditional invariance under scaling requires third order deriva- 
tives. Derivatives are by nature translationally invariant. Ro- 
tational invariance in 2-d is achieved by using rotationally 
symmetric operators. 



Megaplus XRC camera 



Martin ffl reports the settings of 7 = 0.45, 0.50, 0.60 for the Kodak 
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2 The Invariants 

The key idea for the design of the proposed invariants is 
to form suitable ratios of the derivatives of the image func- 
tion such that the parameters describing the transformation 
of interest will cancel out. This idea has been used in [|l2| 
to achieve invariance under linear brightness changes, and it 
can be adjusted to the context of gamma correction by - at 
least conceptually - considering the logarithm of the image 
function. For simplicity, we begin with 1-d image functions. 

2.1 Invariance under Gamma Correction 

Let f(x) be the image function, i.e. the original signal, as- 
sumed to be continuous and differentiable, and f 1 (x) = 
p fix) 1 the corresponding gamma corrected function. Note 
that f(x) is a special case of / 7 (x) where 7 = p = 1. Tak- 
ing the logarithm yields 



=ln(p/(x) 7 ) = mp + 7ln/(x) 



(2) 



with the derivatives fL(x) — 7 f'(x)/f (x), and f!!(x) = 
7 (f(x) f"{x) - f'(x) 2 )/f(x) 2 . We can now define the in- 
variant 8i2 7 under gamma correction to be 



e 127 (/(*)) =fg 
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f(x) /"(*)-/'(*)» 



(3) 



The factor p has been eliminated by taking derivatives, and 
7 has canceled out. Furthermore, 0i2 7 turns out to be com- 
pletely specified in terms of the original image function and 
its derivatives, i.e. the logarithm actually doesn't have to be 
computed. The notation Oi2 7 (/(a0) indicates that the in- 
variant depends on the underlying image function f(x) and 
location x - the invariance holds under gamma correction, 
not under spatial changes of the image function. 

A shortcoming of 8i2 7 is that it is undefined where the 
denominator is zero. Therefore, we modify 0i2 7 to be con- 
tinuous everywhere: 



e 





//' 
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ff"-f' 2 
//' 



if//' = 0A//"-/' 2 =0 
if|//'l < \.ff"-.f' 2 \ 



else 



(4) 



where, for notational convenience, we have dropped the 
variable x. The modification entails —1 < 9 m i2 7 < 1. 
Note that the modification is just a heuristic to deal with 
poles. If all derivatives are zero because the image function 
is constant, then differentials are certainly not the best way 
to represent the function. 



2.2 Invariance under Gamma Correction and 
Scaling 

If scaling is a transformation that has to be considered, then 
another parameter a describing the change of size has to 
be introduced. That is, scaling is modeled here as variable 
substitution the scaled version of f(x) is g(ax) — 

g(u). So we are looking at the function 

/ 7 (x) =ln(p/(a;) 7 ) = Inp + jlng(ax) = g 7 (u) 

where the derivatives with respect to x are g 7 (it) = 
7 ag'(u)/g(u), g>>(u) = ja 2 (g(u) g"(u) - g'(u) 2 )/g(u) 2 , 
and g'^u) = 7 a 3 (g"'(u)/g(u) - 3 g> '(«) g" \u) / 'g(u) 2 + 
2 g' (u) 3 1 g(u) 3 ). Now the invariant Oi23 7 (<7(w)) is obtained 
by defining a suitable ratio of the derivatives such that both 
7 and a cancel out: 



0i23 7 (<?(«)) - g,, {u) 2 



g 2 gig"'-3gg a g"+2 a ' 4 



(5) 



Analogously to eq. (H), we can define a modified invariant 



e 



ml23 7 






if condl 


g 2 g'g"'-3gg' 2 g"+2g' 4 


if cond2 


gZg" i -2gg>^g"+ S ' 4 


g 2 g" 2 -2gg' 2 g"+ g' 4 


else 


g 2 g'g"'-3gg' 2 g"+2g' i 



(6) 



where condition condl is g 2 g' g 1 " — 3 g g' g" + 2 g' = 



A g*g"- - 2gg' 2 g" 



■>....' " 1 g' A = 0, and condition cond2 is 

.2 //2 



\9 2 9' 9"' ••'"/'/ //" • 2'/" < ISV" - 2 gg"g» 
Again, this modification entails — 1 < @ m i23 7 

< 1. 

2.3 An Analytical Example 

It is a straightforward albeit cumbersome exercise to verify 
the invariants from eqs. (^) and (^) with an analytical, dif- 
ferentiable function. As an arbitrary example, we choose 

f(x) — 3xsin(27rx) + 30 

The first three derivatives are f'(x) = 3sin(27ra;) + 
6irx cos(27rx), f"(x) = 127rcos(27r.-r) — \2ir 2 x sin(27ra;), 
and f"'{x) = -367r 2 sin(27ra;) - 247r 3 xcos(27ra;). Then, 
according to eq. 0i2 7 (/(aO) = (3xsin(27ra;) + 
30) (3x sin(27ra;) + 67rx cos(27rx)) / ((3ccsin(27ra:) + 
30) (127rcos(27rx) - 127r 2 xsin(27ra;)) - (3sin(27ra;) + 
6irx cos(27rcc)) 2 ). 

If we now replace f(x) with a gamma corrected version, 
say /o.45(^) = 255 1 - - 45 • 3xsin(27ra;) + 30) - 45 , 
the first derivative becomes /o. 45 0*0 = 255 55 • 
0.45 (3sin(27ra;) + 30)-° 55 (3 sin(27ra;) + 6irx cos(27ra;)), 
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the second derivative is jo'45 (x) = — 255 55 • 
0.45 • 0.55(3sin(27rx) + ' 30)- 155 (3 sin(27rx) + 
67racos(27ra;)) 2 + 255 55 • 0.45 (3a; sin(2vra;) + 
30)-°- 55 (127rcos(27ra:) - 12tt 2 x sin(27nr)), and the third 
is /o.4 5 (X) = 255 - 55 • 0.45 (3 sin(27ra;) + 30)-°- 55 (1.55 • 
0.55 (3 sin(27ra; ) + 30)" 2 (3sin(27nr) + 67racos(27ra;)) 3 + 
(-3) • 0.55(3sin(27rx) + 30)^(3 sin(27rx) + 
67rx cos(27rx)) (127r cos(27rx) — 12-k 2 x sin(27ra;)) + 
(— 367r 2 sin(27ra;) — 24n 3 x cos(2tt2;))). If we plug 
these derivatives into eq. (|J), we obtain an expression for 
@i2 7 (/o.45(x)) which is identical to the one for Oi2 7 (/(a;)) 
above. The algebraically inclined reader is encouraged to 
verify the invariant 0i23 7 for the same function. 




Figure 2: An analytical example function, (left) f(x) — 
3x sin(27nz:) + 30; (right) f 7 (x) = pf(x)f, 7 = 0.45. 
(first row) original functions (second row) first derivatives; 
(third row) second derivatives; (fourth row) third deriva- 
tives; (fifth row) O m i2 7 ; (sixth row) O m i23 7 . 

Fig. U shows the example function and its gamma cor- 
rected counterpart, together with their derivatives and the 
two modified invariants. As expected, the graphs of the in- 
variants are the same on the right as on the left. Note that 



the invariants define a many-to-one mapping. That is, the 
mapping is not information preserving, and it is not possible 
to reconstruct the original image from its invariant represen- 
tation. 

2.4 Extension to 2-d 

If m i2 7 or m i23 7 are to be computed on images, then 
eqs. (|^) to (^) have to be generalized to two dimensions. 
This is to be done in a rotationally invariant way in order 
to achieve invariance under similarity transformations. The 
standard way is to use rotationally symmetric operators. For 
the first derivative, we have the well known gradient magni- 
tude, defined as 

V(x,y) = yjiz + I*:=l' (7) 

where I(x, y) is the 2-d image function, and I x , I y are par- 
tial derivatives along the x-axis and the y-axis. For the sec- 
ond order derivative, we can use the linear Laplacian 

V 2 (x,y) = I xx + I yy :=I" (8) 

Horn [Q] also presents an alternative second order derivative 
operator, the quadratic variation^ 

QV(x, y) = yjl* x + 2Il y + P yy (9) 

Since the QV is not a linear operator and more expensive to 
compute, we use the Laplacian for our implementation. For 
the third order derivative, we can define, in close analogy 
with the quadratic variation, a cubic variation as 

CV(X, y) = \/ I X XX 3Ix Xy + 3I xyy + lyyy - = I 

(10) 

The invariants from eqs. (fy to (^) remain valid in 2-d if 
we replace /' with I', f" with I", and /"' with I'". This can 
be verified by going through the same argument as for the 
1-d functions. Recall that the critical observation in eq. (||) 
was that 7 cancels out, which is the case when all derivatives 
return a factor 7. But such is also the case with the rotation- 
ally symmetric operators mentioned above. For example, if 
we apply the gradient magnitude operator to I(x, y), i.e. to 
the logarithm of a gamma corrected image function, we ob- 
tain 

returning a factor 7, and analogously for V 2 , QV, and CV. 
A similar argument holds for eq. (Q) where we have to show, 
in addition, that the first derivative returns a factor a, the 
second derivative returns a factor a 2 , and the third derivative 
returns a factor a 3 , which is the case for our 2-d operators. 

3 Actually, unlike Horn, we have taken the square root. 



3 



2.5 Differential Operators 



3 Experimental Data and Results 



While the derivatives of continuous, differentiable func- 
tions are uniquely defined, there are many ways to imple- 
ment derivatives for sampled functions. We follow Schmid 
and Mohr JTlj], ter Haar Romeny [13J, and many other re- 
searchers in employing the derivatives of the Gaussian func- 
tion as filters to compute the derivatives of a sampled image 
function via convolution. This way, derivation is combined 
with smoothing. The 2-d zero mean Gaussian is defined as 



G = 



1 



2vro- 2 



(ID 



The partial derivatives up to third order are G x = —xja 2 G, 
G u = -y/a 2 G, G xx = (x 2 - a 2 ) /a 4 G, G xy = xy/a* G, 



(y 2 - a 2 )/a 4 G, G xxx 
(a 2 y-x 2 y)/a 6 G, G xy 



(3a 2 x - x 3 )/a 6 G, 
(o 2 x - xy 2 )/a 6 G, 



G XX y \y Lf a* y J J yj , \-J X yy 

Gyyy = (3ct 2 j/ — y 3 ) / a 6 G. They are shown in fig. |[ We 
used the parameter setting a = 1.0 and kernel size 7x7. 
With these kernels, eq. (ft), for example, is implemented as 



O127 = 



I^/(I®G X ) 2 + {I®Gy) 2 



1(1® G XX + I ® Gyy) 



((I®G X ) 2 + {I®Gy) 2 ) 



at each pixel (x, y), where © denotes convolution. 






Figure 3: Partial derivatives of the 2-d Gaussian, (first row) 
G x , G y , G xx ; (second row) G xy , G yy , G xxx ; (third row) 

G XX y, G X yy, Gyyy. 



We evaluate the invariant 6 m i2 7 from eq. (Q) in two differ- 
ent ways. First, we measure how much the invariant com- 
puted on an image without gamma correction is different 
from the invariant computed on the same image but with 
gamma correction. Theoretical, this difference should be 
zero, but in practice, it is not. Second, we compare tem- 
plate matching accuracy on intensity images, again without 
and with gamma correction, to the accuracy achievable if in- 
stead the invariant representation is used. We also examine 
whether the results can be improved by prefiltering. 

3.1 Absolute and Relative Errors 

A straightforward error measure is the absolute error, 



^Gc(i,j) = |©GC(«, j) ~ ©0GC(«, i)| 



(12) 



where "OGC" refers to the image without gamma correction, 
and GC stands for either "SGC" if the gamma correction is 
done synthetically via eq. ([!]), or for "CGC" if the gamma 
correction is done via the camera hardware. Like the in- 
variant itself, the absolute error is computed at each pixel 
location (i, j) of the image, except for the image boundaries 
where the derivatives and therefore the invariants cannot be 
computed reliably. 




Figure 4: Example image WoBA: (a) no gamma correction, 
"OGC"; (b) gamma correction by camera, "CGC"; (c) syn- 
thetic gamma correction, "SGC". 



Fig. [| shows an example image. The SGC image has 
been computed from the OGC image, with 7 = 0.6. Note 
that the gamma correction is done after the quantization of 
the OGC image, since we don't have access to the OGC im- 
age before quantization. 

Fig. ^ shows the invariant representations of the image 
data from fig. Q and the corresponding absolute errors. Since 
— 1 < O m i2 7 < 1, we have < Agp < 2. The dark points 
in fig. ||, (c) and (e), indicate areas of large errors. We ob- 
serve two error sources: 

• The invariant cannot be computed robustly in homoge- 
neous regions. This is hardly surprising, given that it 
is based on differentials which are by definition only 
sensitive to spatial changes of the signal. 
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20 40 60 SO 100 120 20 40 BO 80 100 120 



(e) (f) 

Figure 5: Absolute errors for invariant m i2 7 , no pre- 
filtering. (a) image WoBA, 0GC; (b) 9 gc; (c) A C gc; 
(d) QcgC, (e) A SGC ; (f) @sgc- 



• There are outliers even in the SGC invariant represen- 
tation, at points of very high contrast edges. They are a 
byproduct of the inherent smoothing when the deriva- 
tives are computed with differentials of the Gaussian. 
Note that the latter put a ceiling on the maximum gra- 
dient magnitude that is computable on 8 -bit images. 

In addition to computing the absolute error, we can also 
compute the relative error, in percent, as 

8cGc(i,j) = 100 A C Gc(h j) / ©ogc(«, j) 

(13) 

Then we can define the set RP e of reliable points, relative 
to some error threshold e, as 

RP e = {(i,j)\6(i,j)<e} (14) 

and PRP e , the percentage of reliable points, as 

PRP e = 100 |RP e jn (15) 




Figure 6: Reliable points RP C for invariant m i2 7 , in black, 
image WoBA, without and with prefiltering. (a) a pre =Q, e= 
5.0; (b)(7 pre =0, e=10.0; (c) a pre =0, e=20.0; (d) a pre =l.0, 
e= 5.0; (e) a pre =l.0, e=10.0; (f) a pre =l.0, e=20.0. 



where n is the number of valid, i.e. non-boundary, pixels in 
the image. Fig. ^| shows, in the first row, the reliable points 
for three different values of the threshold e. The second row 
shows the sets of reliable points for the same thresholds if we 
gently prefilter the 0GC and CGC images. The correspond- 
ing data for the ten test images from fig. |ll] is summarized 
in table 0. 



image 


5.0 


10.0 


20.0 


5.0 


10.0 


20.0 


Build 


13.3 


24.9 


43.8 


16.0 


29.5 


49.3 


WoBA 


15.6 


29.0 


48.2 


19.0 


35.7 


58.9 


WoBB 


16.5 


28.7 


47.1 


21.4 


37.7 


58.1 


WoBC 


18.5 


33.6 


53.5 


24.0 


41.4 


65.3 


WoBD 


13.0 


23.9 


41.9 


16.7 


32.6 


55.6 


Cycl 


15.4 


28.3 


45.9 


22.6 


38.8 


57.6 


Sand 


14.5 


27.2 


44.7 


22.0 


38.5 


57.6 


ToolA 


5.6 


10.7 


20.1 


7.4 


14.7 


27.1 


ToolB 


6.1 


12.0 


22.7 


8.3 


15.7 


28.6 


ToolC 


5.6 


11.1 


20.8 


7.9 


15.1 


28.3 


median 


13.9 


26.1 


44.3 


17.9 


34.2 


56.6 


mean 


12.4 


22.9 


38.9 


16.5 


30.0 


48.6 



Table 1: Percentages of reliable points for O m i2~f, CGC im- 
ages, for e=5.0, 10.0, 20.0. The three numerical columns 
on the left show PRP £ without prefiltering, the three right 
columns with Gaussian prefiltering, a pre =1.0. 

Derivatives are known to be sensitive to noise. Noise can 
be reduced by smoothing the original data before the invari- 
ants are computed. On the other hand, derivatives should 
be computed as locally as possible. With these conflicting 
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goals to be considered, we experiment with gentle prefilter- 
ing, using a Gaussian filter of = 1.0. The size of the 

Gaussian to compute the invariant Q m i2-y is set to c<ier = l-0. 
Note that a pre and ader can not be combined into just one 
Gaussian because of the non-linearity of the invariant. 

With respect to the set of reliable points, we observe that 
after prefiltering, roughly half the points, on average, have 
a relative error of less than 20%. Gentle prefiltering consis- 
tently reduces both absolute and relative errors, but strong 
prefiltering does not. 



3.2 Template Matching 

Template matching is a frequently employed technique in 
computer vision. Here, we will examine how gamma correc- 
tion affects the spatial accuracy of template matching, and 
whether that accuracy can be improved by using the invari- 
ant m i2 7 - An overview of the testbed scenario is given in 



OGC intensity 



CGC intensity 







search 


template 


(correlation) 


i 


0y 














search 


template 


(correlation) 




OGC invariant 



CGC invariant 



Figure 7: The template location problem: A query template 
is cut out from the OGC intensity image and correlated with 
the corresponding CGC intensity image. We test if the cor- 
relation maximum occurs at exactly the same location as in 
the OGC intensity image. The same process is repeated with 
the invariant representations of the OGC and CGC images. 



fig. [7]. A small template of size 6x8, representing the search 
pattern, is taken from a OGC intensity image, i.e. without 
gamma correction. This query template is then correlated 
with the corresponding CGC intensity image, i.e. the same 
scene but with gamma correction switched on. If the cor- 
relation maximum occurs at exactly the location where the 
OGC query template has been cut out, we call this a correct 
maximum correlation position, or CMCP. 

The correlation function s(x,y) employed here is based 
on a normalized mean squared difference c(x, y) [Q]: 

s = max(0, 1 — c) 

£ id ((I(x + i,y + j)-I)-(T(i,j)-T)f 
c = 

\IT~Mx + i,y + j)- I) 2 LMij - T? 



where I is an image, T is a template positioned at (x,y), 
I(x, y) is the mean of the subimage of / at {x, y) of the same 
size as T, T is the mean of the template, and < s < 1. The 
template location problem then is to perform this correlation 
for the whole image and to determine whether the position 
of the correlation maximum occurs precisely at {x, y). 









m. ^ 

a- 




! J 




Figure 8: Matched templates, image WoBA: (left) intensity 
data; (right) invariant representation. Black box=query tem- 
plate, white box=matched template. 




Figure 9: Correlation matrices, image WoBA. (left) intensity 
data; (right) invariant representation. 

Fig. ^| demonstrates the template location problem, on 
the left for an intensity image, and on the right for its in- 
variant representation. The black box marks the position of 
the original template at (40,15), and the white box marks 
the position of the matched template, which is incorrectly 
located at (50,64) in the intensity image. On the right, the 
matched template (white) has overwritten the original tem- 
plate (black) at the same, correctly identified position. Fig. ^ 
visualizes the correlation function over the whole image. 
The white areas are regions of high correlation. 

The example from figs. [| and ^| deals with only one arbi- 
trarily selected template. In order to systematically analyze 
the template location problem, we repeat the correlation pro- 
cess for all possible template locations. Then we can define 
the correlation accuracy CA as the percentage of correctly 
located templates, 



CA, 



= 100 I CMCP, 



(16) 



where tn x tm is the size of the template, CMCPt„ xtm 
is the set of correct maximum correlation positions, and n, 
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1 m Q CT£1 

Ullage 


Tnt/fl 
1I1L/U 


Tnt/l n 


Tnv/f) 
1I1V/U 


Tnv/1 n 
1I1V/ 1 .u 


Dui x a 


OJ.U 


/ O.U 


OJ.O 


07.J 


Mara 


55.5 


45.0 


75.7 


80.4 


WO .Do 


J7.J 




S9 7 


^7 fi 


WoBC 


67.2 


58.3 


68.9 


78.7 


WoBD 


31.6 


29.2 


48.0 


67.4 


Cycl 


60.5 


45.4 


98.6 


99.4 


Sand 


50.5 


40.9 


85.2 


94.4 


ToolA 


41.7 


35.3 


60.2 


68.0 


ToolB 


29.5 


23.4 


45.7 


54.1 


ToolC 


42.1 


27.8 


42.5 


48.4 


median 


46.3 


38.1 


64.6 


73.4 


mean 


50.3 


41.4 


66.3 


73.8 



Table 2: Correlation accuracies CA, template size 6x8, left 
columns for intensity data, right columns for the invariant 
representation, without and with prefiltering. 



Figure 10: Binary correlation accuracy matrices, white 
pixels=CMCP6x8> image WoBA. (a) intensity image, 
a pre =Q; (b) intensity image, a pre = 1.0; (c) invariant rep- 
resentation, <j pre =0; (d) invariant representation, a pre =l.Q. 



again, is the number of valid pixels. We compute the cor- 
relation accuracy both for unfiltered images and for gently 
prefiltered images, with a pre = 1.0. Fig. |To| shows the bi- 
nary correlation accuracy matrices for our example image. 
The CMCP set is shown in white, its complement and the 
boundaries in black. We observe a higher correlation accu- 
racy for the invariant representation, which is improved by 
the prefiltering. 
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Figure 11: Test images: (a) Build; (b) WoBA; (c) WoBB; 
(d) WoBC; (e) WoBD; (f) Cycl; (g) Sand; (h) ToolA; 
(i) ToolB; (j) ToolC. 



Figure 12: Correlation accuracies, (top) template size 6x8; 
(bottom) template size 10 x 10. The entries at x=l refer to 
Build, at x=2 to WoBA, etc. (circles, lower line) intensity 
images; (stars, center line) invariant representation, a pr e=0; 
(diamonds, upper line) invariant representation, <j pre =l.0. 
The markers on the left hand side indicate the means, the 
markers on the right hand side the medians. 
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We have computed the correlation accuracy for all the 
images given in fig. [TT] . The results are shown in table ^ and 
visualized in fig. |12| We observe the following: 

• The correlation accuracy CA is higher on the invariant 
representation than on the intensity images. 

• The correlation accuracy is higher on the invariant rep- 
resentation with gentle prefiltering, a pre = 1.0, than 
without prefiltering. We also observed a decrease in 
correlation accuracy if we increase the prefiltering well 
beyond a pre = 1.0. By contrast, prefiltering seems to 
be always detrimental to the intensity images CA. 

• The correlation accuracy shows a wide variation, 
roughly in the range 30%. . . 90% for the unfiltered in- 
tensity images and 50%. . . 100% for prefiltered invari- 
ant representations. Similarly, the gain in correlation 
accuracy ranges from close to zero up to 45%. For our 
test images, it turns out that the invariant representation 
is always superior, but that doesn't necessarily have to 
be the case. 

• The medians and means of the CAs over all test images 
confirm the gain in correlation accuracy for the invari- 
ant representation. 

• The larger the template size, the higher the correlation 
accuracy, independent of the representation. A larger 
template size means more structure, and more discrim- 
inatory power. 

4 Conclusion 

We have proposed novel invariants that combine invariance 
under gamma correction with invariance under geometric 
transformations. In a general sense, the invariants can be 
seen as trading off derivatives for a power law parameter, 
which makes them interesting for applications beyond im- 
age processing. The error analysis of our implementation 
on real images has shown that, for sampled data, the invari- 
ants cannot be computed robustly everywhere. Neverthe- 
less, the template matching application scenario has demon- 
strated that a performance gain is achievable by using the 
proposed invariant. 
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